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ABSTRACT 

An  analysis  was  performed  to  determine  range  acc- 
uracy for  echo  ranging,  bottom-bounce  sonar  systems.  The 
model  for  analysis  assumes  isovelocity  conditions  and  a 
sloping  bottom.  The  range  to  the  target  is  determined  from 
a knowledge  of  the  echo  travel  time  and  measurement  of  the 
required  environmental  parameter  (e.g.,  water  depth,  bottom 
slope).  The  results  indicate  that  the  range  error  due  to 
bottom  slope  is  independent  of  range  for  small  bottom  slopes. 
The  range  errors  due  to  inaccuracies  in  measuring  the  bottom 
depth  and  slant  range  are  also  examined  and  the  rms  error 
in  range  is  determined  for  a representative  environment. 
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SECTION  I 
INTRODUCTION 

For  effective  utilization  of  a sonar  system,  the 
user  must  estimate  the  target  range  from  a knowledge  of  the 
environmental  parameters  and  measurement  of  the  echo  travel 
time.  In  addition,  an  estimate  of  the  range  accuracy  is  re- 
quired if  the  sonar  detection  is  to  be  used  as  a basis  for  a 
fire  control  solution^. 

In  this  paper  an  analytic  expression  is  developed 
for  the  determination  of  target  range  in  an  isovelocity  med- 
ium with  a sloping  bottom,  given  the  depth  of  the  water,  the 
slope  of  the  bottom  and  the  travel  time  for  the  echo.  Ex- 
pressions are  subsequently  developed  for  the  range  errors  re- 
sulting from  the  uncertainties  in  the  measurement  of  these 
parameters  and  representative  RMS  range  errors  are  calcul- 
ated for  various  ranges  and  pulse  lengths. 

(2) 

The  work  by  T.N.  Reynolds'  1 is  extended  by  con- 
sidering errors  in  sound  velocity,  bottom  depth  determina- 
tion, and  transmission  time  measurement.  In  this  paper  a 
two-dimensional  model  is  considered,  and  the  expression  ob- 
tained herein  reduces  to  Reynolds'  result  for  small  bottom 
slope.  The  rate  of  change  of  the  range  with  depth  and 
slant  range  (travel  time)  are  also  evaluated. 

Figure  1 illustrates  the  geometry  for  the  analysis. 
In  the  figure: 

L = slant  range  to  target  and  equals  cT/2,  where 
c is  the  velocity  of  sound  and  T is  the  travel 
time, 

d = perpendicular  distance  from  sonar  to  the  bottom 
of  the  ocean, 

d'=  vertical  depth  to  the  bottom  at  point  of  re- 
flection. 


1. 
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0 s bearing  angle  measured  from  the  sea  surface, 
fx=  angle  of  bottom  slope,  positive  as  shown, 

R = horizontal  range  from  sonar  to  target, 

distance  from  sonar  to  point  of  specular  re- 
flection on  bottom, 

distance  from  specular  reflection  on  bottom  to 
target, 

Xj=  horizontal  distance  from  sonar  to  specular  re- 
flection point,  and 

x^  horizontal  distance  from  specular  reflection 
point  to  target. 

The  difference  between  the  depth  of  the  sonar  and 
the  depth  of  the  target  is  small  compared  to  the  total  depth 
of  the  ocean.  Consequently,  the  error  introduced  by  assuming 
that  both  are  at  the  surface,  for  deep  bottoms,  is  negligible. 
In  the  error  analysis  later  on,  their  finite  depths  will  be 
considered  in  greater  detail. 
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SECTION  II 

DEVELOPMENT  OF  THE  RANGE  FORMULA 


The  slant  range  L is  given  by 

L - + i2  ( 

and  the  horizontal  target  range  is  seen  to  be 

R - + x2  ( 

By  the  Pythagorean  Theorem 

<V2  - xi2  + d'2  ( 

(i2)2  = x22  + d'2  ( 

Referring  to  Figure  1,  the  following  expressions  are  easily 
written: 


sin  0 

cos  0 


*1 

xi 

*1 


and. 


sin  (0+a)  = -s— 
*1 

sin  (0+a)  = 4— 
l2 


Subtracting  (4)  from  (3)  gives 

V ~ v*2 


(^)2  - U9)2  - (*i)2 


( 

( 

( 

( 

0 


Now  substituting  (5)  and  (8)  into  (1)  yields 
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Equations  (5)  and  (7)  may  be  solved  for  d' 

, i _ d sin  9 
sin  (0+a) 

Substituting  (11)  into  (10)  gives 

t _ d sin  9 f 1 , 1 1 

sin (04 a)  Lsin0  sin  (9+2  a)  J 

utilization  of  the  identity:  sin  A + sin  B = 2 [sin 

results  in 


(11) 


(12) 


sin  (0+2 a)  = -j-  cos  a 


(13) 


By  the  Law  of  Sines,  an  expression  for  R in  terms  of  L,  0,  and 
a may  be  obtained  from  Figure  1 as  follows: 


^2  _ R 

sin  9 = sin (20+ 2 a) 


(14) 


sin(04-2a) 


R 

sin  (2W+Za) 


(15) 


or. 


- R sin  9 

^2  = sin(20+2a) 


R sin 


H = sin(2£+2^T 


(16) 

(17) 


Sin  L = + A 2 one  obtains 

L - sln(mST  (sln  9 + sln  <9+2°>]  (l8) 

but,  sin  0 + sin  (0+2a)  = 2 sin  (0+a)  cos  a by  the  expansion 
for  the  sin  (A+B),  and  sin  (20+2a)  * 2 sin  (0+a)  cos  (0+a)» 


5. 
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combining 


a.  L 

b. 


R COS  q 

cos(9+a) 


R _ cos  JHsd 
L cos  a 


(19) 


Now  expanding  the  left  side  of  equation  (13), 

sin  j[(0-ta)]  + sin  (9+a)  cos  a + cos  (9+o)  sin  a (20) 

Divide  both  sides  of  equation  (13)  by  cos  a and  substitute  eq- 
uation (19b),  then 


sin(9+a)  + ^ sin  as-^  (21) 

Solving  for  sin(9+a)  yields 

sin (9+o ) = Vl  - cos5(0+a)  = ~ - £ sin  a (22) 

o 2 

But,  from  (19b)  — 7 cos  a = cos  (9+a) 

L 

Hence,  substituting  into  (22),  squaring  both  sides  and  col- 

2 

lecting  terms,  yields  an  expression  in  R 


R2  - 4dR  sino  + (4d2-L2)  = 0 (23) 

Using  the  quadratic  formula  to  solve  for  R and  choosing  the 
positive  sign  of  the  square  root  on  physical  grounds,  results 
in 


R ■ 2d  sin  a +*V4d^sin2o  + (L2-4d^) 


(24) 
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FIGURE  2o  GEOMETRY  FOR  CALCULATION  OF  RANGE 


where : 


R 

L 

2d 


range 

slant  range 

which  is  twice  the  perpendicular  distance  to 
the  bottom  measured  from  the  sonar. 


one  obtains 


and 


a)  sin  0O  = ^7 


b)  cos  0 


L2  -4d2 


(25) 


Hence,  (24)  may  be  written 


R = 2d  sina  + L cos  0 


VT  + tan2  0 sin2  a (26) 


Now,  since  sin  a is  small,  the  binomial  expansion  is  applicable. 
Therefore, 

,n2 

T 


Vl  + tan2  0 sin2  0 ~ 1 + S3,%  - tan2  0, 
o 


But 


si.n_QL  tan2  9 « -1 

? o 


l 


& 

I 
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In  this  approximation,  (26)  becomes 


or. 


R = 2d  sin  a + L cos  0 

o 


R L cos  0Q  + 2da 


(27) 


since  a is  small,  sin  a ^ a. 

Expression  (27)  is  composed  of  two  terms:  Rq  = L cos  0Q,  and 
AR  “ 2d  sin  a«  The  RQ  term  is  the  range  for  zero  bottom  slope 
and  the  AR  term  arises  due  to  the  non-zero  bottom  slope.  The 
interesting  observation  is  that  this  term,  interpreted  as  an 
error  (or  correction  factor)  is  independent  of  the  range! 

Figure  3 is  a plot  of  the  range  correction  term, 

AR  = R-R0  vs  bottom  slope  for  a water  depth  of  5 kyd  and  a 
30  kyd  slant  range  to  the  target.  Note  that  the  graph  is  for 
positive  a (depth  decreasing  from  sonar  to  target);  the  range 
correction  term  becomes  negative  for  negative  a« 


8. 
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SECTION  III 

DETERMINATION  OF  RANGE  DERIVATIVES 

From  equations  (13)  and  (19)  the  range  changes  with 
respect  to  a,  L and  d may  be  found. 

Since  R = R(L,  d,  a),  using  the  chain  rule. 


dR  = Sr  + aR  ae 

da  oa  da  "5a 


3£ 


and 


sin(0+2«) 

2d 

0 

eina 

+ 2m 

sin(0+a)  1 

cos  q 

L 

cos(0+2a) 

ji 

cosa  | 

dR  . 

. dR 

dR  c 

)0 

HE  “ 

■ 31 

+ d§  c 

5E 

ii 

cos(9+a) 

cosa 

+ ¥ 

sin 

cos 

dR 

dR 

. dR 

de 

HH  ‘ 

■ d 3 

33 

dR 

HH  - 

- 2 

sin(0+a) 
cos  (d'+Za) 

(29) 


(30) 


(Where  it  is  understood  that  the  "total"  derivatives  are  really 
partial  derivatives  with  the  appropriate  parameters  held  con- 
stant) . 

For  a = 0 equations  (28)  to  (30)  reduce  to: 

•gS  = L sin  ©q  = 2d  (The  rate  of  change  of  (31) 

R with  respect  to  a is 
independent  of  R for 
fixed  travel  time!) 


( 


10. 


- sec  e (32) 

H3  = -2  tan  9 (33) 


The  result  of  equation  (31)  is  contrary  to  intuition  but  may  be 
more  readily  seen  by  inspecting  equation  (27). 
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SECTION  IV 

CALCULATION  OF  RMS  ERROR  IN  RANGE 
Consider  the  error  due  to  errors  in  a,  L,  and  d. 

AR  = fr  + + H ad  <*» 

Since  the  differentials  are  independent,  the  expressions  for 
the  RMS  error  in  range  is 

oR  =V<ARa)2  + (arl)2  + (ARd)2  (35) 

To  obtain  some  representative  estimates  of  range  errors  which 
could  occur,  calculations  were  performed  for  the  environment 
values  listed  in  Table  1. 


TABLE  I 


Table  of  Assumed  Environmental  Values 


Quantity 

Value 

Uncertainty 

Speed  of  Sound(fps) 

5000 

2.5(.05%)(3,4) 

Water  Depth  (yd) 

5000 

25(5) 

Target  depth  (yd) 

- 

150 

Bottom  Slope  (deg) 

0 

X-S 

CM 

] 

• 

rH 

12. 
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SECTION  V 

SLANT  RANGE  ERRORS 


0 


! 


There  are  two  sources  of  error  attributable  to  the 
velocity  profile.  The  first  Is  due  to  limitations  In  the  pre- 
cision of  measurement  apparatus,  and  the  second  Is  due  to  hori- 
zontal velocity  gradients  (i.e.,  the  velocity  profile  varies 
with  range  and  bearing).  The  exact  magnitude  of  the  range 
errors  associated  with  these  parameters  uncertainties  is  not 
known.  In  the  isovelocity  model,  it  is  felt  that  this  source 
of  error  may  be  adequately  represented  as  an  uncertainty  in 
the  mean  sound  velocity. 

Since 

t cT 

L " T 


(38) 


dL  = dt  + i dc 


and  the  error  in  L,  AL  is  given  by 


AT  * ^ AT  4-  AC 

AL  = 7 AT  + y y 


(39) 


or , equiva len t ly , 


AL  = L y + y 


(40) 


(41) 


where  the  uncertainty  in  travel  time  is  simply  the  pulse  length, 
t,  for  CW  signals  (or,  for  linear  FM  signals,  the  resolution 
is  given  by  the  reciprocal  of  the  bandwidth) . The  uncertainty 
in  average  velocity  is  estimated  to  be  .05%  (y  = .0005). 


For  R = 15  Kyds.  and  t **  10  msec,  one  obtains  (see 


Figure  2): 


■Vr2  + 


4D 


=V(15)2  + 4(5)2  = 18.02, 


(42) 


13 
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and 


- sin 


-1  10 


15707- 


- 33°  40 


From  equation  (32),  we  get 


arL  * 


AL  sec  9 

ft 


sec  9 


dc  c t1 

- + tJ 

.02  x 10^5  x 10'f  +£ 


x 10' 


14.6  yds. 


(44) 


1/2 


sec  33 


r v 
[ 


14 
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SECTION  VI 
BOTTOM  DEPTH  ERRORS 

To  estimate  the  range  error  associated  with  depth 
uncertainties,  one  must  consider: 

(1)  the  depth  to  the  bottom  of  the  ocean, 

(2)  the  depth  of  the  sonar,  and 

(3)  the  depth  of  the  target. 

Note  that  in  this  study  the  water  depth  is  measured  perpendicu- 
lar to  the  bottom  (see  Figure  1) . Figure  4 is  similar  to 
Figure  1 with  the  addition  of  finite  sonar  and  target  depths. 

The  virtual  depth  of  the  target  is  now  seen  to  be: 

2d  - (d  + d ) sec  a. 

W S 

For  convenience,  define 


D - d - 1/2  (d  + d ) sec  a 

W S 

To  find  the  errors  in  the  depths,  take  the  derivative  of  D* 


dD  - dd  - 1/2  (ddfc+  dd^)  sec  a + 1/2  (dt  + dg)  tan  a sec  a (46) 

Since  for  most  cases  of  interest  the  sonar  depth  error 
is  negligible,  the  dd  term  may  be  dropped.  For  a zero  mean 

O 

bottom  slope,  equation  (46)  reduces  to 

dD*  - dd  - 1/2  ddt  (47) 

and  since  the  errors  are  independent, 

AD*  -V(dd)2  + (-l/2ddt)2  (48) 

Now  from  equations  (33)  and  (48)  the  error  in  range  due  to 


15. 
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0 


depth  errors  is  found  to  be: 

ARd*  - - 2AD*  tan  0 (49) 

Substituting  the  values  from  the  table 


ARp  - -2V(25)i  + (1/2  150) Z tan  33°  40' 
» 105  yds. 
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SECTION  VII 

ERROR  IN  RANGE  DUE  TO  BOTTOM  SLOPE 

A mean  bottom  slope  of  1-1/2°  is  assumed  based  on 
Reynolds' measurements.  This  is  probably  pessimistic, 
since  he  apparently  attributed  the  entire  range  error  which 
he  measured  to  bottom  slope  and  neglected  other  sources  of 
error  as  discussed  in  this  paper. 

From  equation  (31),  the  error  in  range  due  to  bottom 
slope  may  be  written  as: 

ARa  - 2dAa  (50) 

Hence,  substituting  the  values  from  Table  1,  one  obtains 

arq  = fitr (5000)  (1°5)  = 261  yds 
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SECTION  VIII 
RMS  ERROR  IN  RANGE 

Combining  the  errors  found  in  (44),  (49)  and  (50)  and 
substituting  into  (35)  we  get 


oR  = (261)2  + (105)2  + (14. 6) 2 

= 281  yds. 


(51) 


Table  2 gives  the  components  namely,  AR  , ARn  and  AR. 
that  make  up  the  RMS  error  in  range  at  three  different  ranges 
and  for  six  pulse  lengths.  Figure  5 is  a plot  of  range  error 
(or)  vs  pulse  length  (t)  for  constant  target  range,  and  Figure 
6 is  a plot  of  range  error  (aR)  vs  target  range  (R)  for  fixed 
pulse  lengths. 

From  these  figures  one  may  observe  that  for  the  10ms 
and  100  ms  pulse  lengths  there  is  almost  no  difference  in  range 
accuracy  because  the  dominant  source  of  error  is  due  to  the 
bottom  slope  term.  For  the  500  ms  pulse  length,  the  uncertainty 
due  to  the  pulse  length  dominates  the  other  terms  and  causes  a 
significant  loss  in  range  accuracy. 


19. 


FIGURE  5.  RANGE  ERROR  (op)  VS.  PULSE  LENGTH 


TABLE  2 


Pulse  Length 

R 

d 

L 

0 

ARa 

ARd 

A*L 

mm 

10 

ms 

15 

5 

18.02 

33°40 ' 

261 

105 

14.6 

281 

100 

ms 

15 

5 

18.02 

33° 40 ' 

261 

105 

98.6 

298 

200 

ms 

15 

5 

18.02 

33°40  1 

261 

105 

196.4 

343 

300 

ms 

15 

5 

18.02 

33°40 1 

261 

105 

295 

408 

400 

ms 

15 

5 

18.02 

33°40  ' 

261 

105 

392 

482 

500 

ms 

15 

5 

18.02 

33040  • 

261 

105 

490 

565 

10 

ms 

30 

5 

31.62 

18°15  ’ 

261 

52 

18.7 

267 

100 

ms 

30 

5 

31.62 

18° 15  • 

261 

52 

87.4 

280 

200 

ms 

30 

5 

31.62 

18°15  ' 

261 

52 

173.1 

317 

300 

ms 

30 

5 

31.62 

18°15  ' 

261 

52 

259 

371 

400 

ms 

30 

5 

31.62 

18°15  1 

261 

52 

344 

435 

500 

ms 

30 

5 

31.62 

18° 15  ' 

261 

52 

430 

505 

10 

ms 

45 

5 

46.1 

12°32  1 

261 

35 

25.5 

264 

100 

ms 

45 

5 

46.1 

12°32  ' 

261 

35 

87.0 

277 

200 

ms 

45 

5 

46.1 

12°32  ' 

261 

35 

168.9 

313 

300 

ms 

45 

5 

46.1 

12°32  ' 

261 

35 

251 

364 

400 

ms 

45 

5 

46.1 

12°32  ' 

261 

35 

334 

425 

500 

ms 

45 

5 

46.1 

12°32  ' 

261 

35 

419 

495 
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